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Abstract 

We study translational and spin-isospin symmetry restoration for the 
two-flavor chiral quark-loop soliton. Instead of a static soliton at rest we 
consider a boosted and rotating hedgehog soliton. Corrected classical me- 
son fields are obtained by minimizing a corrected energy functional which 
has been derived by semi-classical methods ('variation after projection'). 

We evaluate corrected meson fields in the region 300 MeV < M < 
600 MeV of constituent quark masses M and compare them with the uncor- 
rected fields. We study the effect of the corrections on various expectation 
values of nuclear observables such as the root-mean square radius, the 
axial- vector coupling constant, magnetic moments and the delta- nucleon 
mass splitting. 



1 Introduction 



Chiral soliton models have turned out to be a fruitful approach to the structure of 
light baryons. We consider the soliton of the two-flavor chiral quark-loop model, 
which is equivalent to the semi-bosonized Nambu & Jona-Lasinio (NJL) model, 
with hedgehog meson fields restricted to the chiral circle (for a detailed review 
cf. ref. [0, 0])- Mesons are treated as mean fields (zero boson loop) whereas the 
polarization of the Dirac sea is fully taken into account. In order to obtain a 
system with baryon number B=l one has to add N c =3 valence quarks occupying 
the lowest single-particle level with positive energy. Static, self-consistent mean- 
field solitons are obtained by minimizing the effective energy. For a detailed 
description of procedure, formalism and notation we refer the reader to ref. || . 

As it is well known the mean-field hedgehog soliton breaks both the trans- 
lational and rotational symmetry. The soliton is not an eigenstate of the total 
linear (P) and angular (J) momentum and of the isospin operator (T). As a 
result, the expectation value (^P 2 ^ does not vanish, and i^J 2 ^ and (T 2 ^ do not 
have the values appropriate for a nucleon or a A isobar. The energy functional, 
which defines the solitonic field configuration, is contaminated by spurious con- 
tributions. There are various method to exclude the spurious energies and to de- 
termine the energy of a state with definite quantum numbers. In non-relativistic 
many-particle physics the spurious energies can be obtained within a certain ap- 
proximation using Peierls-Yoccoz Q or Peierls-Thouless || projection techniques 
for the linear and angular momenta. Unfortunately those methods require the 
definition of a Fock state which is very involved if dealing with the polarized 
and regularized Dirac sea. In the case of the Skyrme model it has turned out 
that the consideration of RPA fluctuations around the mean-field solution allows 
a treatment of the corresponding eigenmodes [|, [F| analogously to the case of 
two-dimensional soliton models ||. Recently this approach has also been used 
to calculate the quantum corrections to the effective soliton energy in the chiral 
quark- loop model within some approximations ||. 

In the present paper we follow a simplified approach and rely on the semi- 
classical pushing and cranking approximations ||10|| . Hereby one is considering a 
soliton with the appropriate quantum numbers by pushing the static hedgehog 
with the velocity V and rotating its isospin coordinates with the angular velocity 
Q where V and Q are regarded as Lagrange multipliers. Due to the hedgehog 
symmetry the angular momentum of the soliton is fixed to J=— T and one has 
automatically ^J 2 ^ = (t 2 ^. All the studies which have been done so far per- 
formed the projection after the variation, which means that first the mean field 
configurations are calculated by minimizing the static soliton energy and after- 
wards the rotational and translational correction terms are added. It is the aim 
of this paper to investigate if and how the properties of the chiral quark-loop 
model change if the more accurate procedure of a variation after projection is 
carried out. In this case one has to minimize a soliton energy which includes the 
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correction term. The resulting corrected meson fields respond to the collective 
translation and rotation of the soliton and deviate from the static ones. The 
inertial parameters correspond to those values which minimize the corrected en- 
ergy functional (self- consistent cranking). Expectation values of quark and meson 
observables are also modified. 

In sect. 2, we define the static hedgehog soliton and introduce a corrected 
energy functional for a soliton with definite momentum and isospin. The field 
configuration which minimizes the corrected energy functional is determined in 
sect. 3. Particular attention is drawn to the asymptotic behavior of the meson 
profile. In sect. 4 we investigate the effect of the corrections on solitonic expec- 
tation values. Conclusions are drawn in sect. 5. 

2 Static mean-field and energy corrections 

In mean-field approximation, the SU(2) NJL lagrangian describes u and d quarks 
interacting with classical meson fields a and jr. Restricted to time-independent 
spherical hedgehog configurations ur(i, r) = c(r), fr(t, r) = rr(r)rj and to the 

chiral circle (cr(r) 2 + 7r(r) 2 = M 2 ^j with unit vector r = r/r, r = |f|, the meson 
fields are uniquely determined by the constituent quark mass M and the profile 
function 0(r) =arctan ^Jy. The mean-field energy 

£ mf [e] = £ m [0] + £ q [0] (i) 

is a functional of G and consists of a purely mesonic part 

E m [9] = ^^4tt / r 2 dr [1 - cos 0(r)] , (2) 

Gr J 

and of the quark energy 

E*[Q] = - lim -^TrLn[d T + h]+ lim -Ur Ln [d T + h ) + N c e val 9(£ val ) (3) 

with the single-particle quark hamiltonians 

h = -ia-V + pMe^ e{r) and h = -ia- V + /3M . (4) 

Here, a = f3j, 7= (7 1 ,7 2 ,7 3 ), 75 = i7°7 1 7 2 7 3 7 4 and (3 = 7 are Dirac matrices, 
m is an average current quark mass m = (m u + m&)/2. The Euclidean time- 
coordinate is denoted by r, while f denotes the vector of Pauli matrices. The 
last term in eq. (^|) describes the energy of the valence quarks which have been 
added to ensure a baryon number B = l for the soliton. It vanishes if the meson 
field is so strong that the energy e v3i \ of the valence level is negative. In this 
case the Dirac sea has already B = 1. The trace Tr in eq. (0), which includes 
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functional trace with anti-periodic boundary conditions in the Euclidean time 
interval (— ^y, ^) and matrix trace over Dirac, flavor and color indices, can be 
expressed by a regularized sum over the eigenvalues of the hamiltonians h and 
h . For details cf.ref. 0. 

The classical meson profiles have to minimize the mean-field energy ([I]) 

5Emi (5) 



50 fr) 



leading to the equation of motion 



P(r) 

O(r) = arctan — , (6) 

S(r) 

which is an implicit equation since scalar and pseudoscalar quark densities S(r) 
and P(r) (cf. ref. [£§) are functionals of the meson fields. At large separations 
from the center of the soliton the meson fields approach to their vacuum values 
cr(r — > oo) = (Xo —M and 7r(r ^oo) = 0, while the difference O(0) — ©(oo) is related 
to the baryon number of the soliton. In order to obtain a soliton with baryon 
number B = 1 we assume 

Pt( \ I ~ 1X ^ or r ~ > ® (7\ 
I for r — > oo 

The implicit equation of motion (||) is usually solved iteratively. Starting from 
a reasonable profile 0° one diagonalizes the hamiltonians (§) and gets a set of 
eigenvalues and eigenfunctions which allow to determine the densities S and P. 
The latter are used on the right side of eq. (|[) in order to get an improved profile. 
This procedure is repeated unless the profile remains unchanged within a certain 
accuracy. It has been shown |TT| that the iteration converges for constituent quark 



masses M > 330 MeV. Mesonic fields which fulfill the equation of motion (||) with 
quark fields which diagonalize the corresponding quark hamiltonian @ are called 
self- consistent solutions of the NJL model in mean-field approximation. They 
have to be distinguished from a parameterized meson field with a predetermined 
shape which approximates the self-consistent field and does not fulfill the equation 
of motion (0) in general. 

The soliton as a mean-field solution with hedgehog shape is not an eigenstate 
of total momentum P and isospin operator T. As a result, the expectation values 
P 2 ) and (T 2 ) do not have the values 



P) = and (T 2 )=T{T+1) (8) 

leon and T— | for a A isobar at rest. The mean field dis- 
tinguishes a definite position (the center) of the soliton and the hedgehog defines 
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a privileged direction (t z is maximal in z direction) what obviously contradicts 
translational and rotational symmetry. The mean-field hedgehog soliton performs 
spurious motions (oscillations around the center of mass and in isospace). These 
spurious motions are not only responsible for the deviation of (P 2 ^ and (T 2 ^ 
from the physical values (H) but contribute to most of the expectation values of 
the soliton. 

The energy of the spurious modes can be approximately calculated in the 
framework of the pushing and cranking approaches often used in non-relativistic 
many-particle physics [|K|] . After semi-classical quantization of the rotational 
degrees of freedom the energy of a soliton with the expectation values (|8|) is given 

I'.v m 

-^corr = -^mf — -^trans ~~ ^rot + -^crank (9) 

with the translational zero-point (center-of-mass) energy 



P 



2\ 



^trans 2M 

and its iso-rotational equivalent 



771U . _ \ / hh (](]) 
trans 2M inert 



;^' 2 ) 9 

E Iot - - gj • (11) 

The parameters M mert and / are the Thouless-Valatin inertial parameters of 



translational and rotational motion of the soliton as a whole. As shown in [[13 
the inertial mass M mert coincides with the static mean- field energy (|Ij). The 
expectation value ^P 2 ^ of the mean-field hedgehog soliton consists of contri- 
butions from valence and sea quarks and is a functional of the meson profile. It 
diverges and therefore has to be regularized. We apply Schwinger's proper-time 



scheme [|4] and relate the cut-off parameter to the constituent quark mass by 



means of the weak pion decay. For details see refs. |0, [15|, [16], [T 



The expectation value of (T 2 ^ is finite. In contrast to the corresponding 
value for the linear momentum it is independent of the shape of the profile func- 
tion and consists of the single-particle expectation values of the 3 valence quarks 
(P 2 ') hh = 3 • \ (| + l) = f • The Dirac sea does not contribute to the isospin of the 
hedgehog soliton. The iso-rotational moment of inertia I is given by a regularized 
Inglis formula and has to be calculated numerically. A term like ( ]TT| ) has a 
well known analog in nuclear physics, namely the band-head energy which has to 
be subtracted from the mean field energy if the cranking procedure is applied to 
a nucleus with spin different from zero. The last term in eq. (^|) 

F* T{T+l) , . 

-^crank ~ 21 
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is the genuine cranking term and results from the global quantized rotation which 
was introduced to obtain the correct value for the isospin. Both the rotational 
corrections (|ll|) and fll2|) can be considered as the result of a single cranking pro- 
cedure in the course of which the expectation value (T 2 ^ has been changed from 
| (resulting from 3 valence quarks) to the physical value The corresponding 
semi-classically quantized angular velocity is given by 



n 



T(T + 1) - (T 



(13) 



where the second term in the numerator stems from the rotational zero-point 
energy of the hedgehog while the factor — i appears because we are working in 
Euclidean space-time. It is worthwhile to mention that the same correction terms 
fllTD and ( |l"2D can be obtained |10[ within an approximate projection method using 
the Kamlah expansion The various contributions (TU-|P2D to the corrected 
soliton energy flj^) are of different order in N c , namely E m { ~ N c , £^ ransrot ~ 1 
and -E^ank ~ even though it will turn out that the cranking term can be 

numerically of the same order of magnitude as the zero-point energies. 



600 



400 



200 



> 

CD 



LU 



-200 



-400 



-600 



-800 




300 



400 



500 
M (MeV) 



600 



Fig. 1: Total energy corrections 

L\Ej — -G/ trans ^rot + ^crank 

for nucleon and A isobar {lower 
part) and their components (up- 
per part) calculated for the static 
soliton profile (broken lines), 
and for the corrected nucleon 
(N) and A profiles (full lines) 
as a function of the constituent 
quark mass M. The cranking 
energy -E^ank °f ^ ne ^ isobar is 
five times larger than the cor- 
responding energy -E^ank of the 
nucleon. The energies of the 
translational and rotational fluc- 
tuations calculated in reference 
H are indicated by O and •, re- 
spectively. 



Fig. 1 illustrates the size of the energy corrections. Each of the terms increases 
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approximately linearly with the constituent quark mass M. The increase of the 
corrections is related to the spatial size of the soliton. At larger constituent masses 
the quarks are stronger bound and the soliton is smaller (cf. fig. 7). Moment of 
inertia and expectation value (-P 2 ) hh behave accordingly. While the moment of 
inertia decreases the momentum fluctuations increase if the size of the soliton is 
reduced at larger M. Already at small quark masses the corrections amount to a 
considerable part of the total mean-field soliton energy of roughly 1230 MeV. At 
M > 480 MeV the total energy correction for the nucleon exceeds one half of the 
soliton energy. For the A isobar the total energy correction is smaller since the 
zero-point corrections, which have to be subtracted, are partly compensated by 
the big positive cranking term. Additionally we have indicated the energies of the 
corresponding mesonic quantum fluctuations calculated in RPA [BJ. Comparing 
them with the semi-classical zero-point energies (|i~0"D and (|ll]) one has to take 
into account that only bound fluctuation have been considered. According to the 
estimate of ref. unbound fluctuations, which appear in the non-confining NJL 
model, contribute roughly 50 percent to the translational modes, while unbound 
rotational modes are negligible. 

The considerable size of the energy corrections in comparison to the mean- 
field energy E m { conflicts with their perturbative calculation. To improve the 
approach we shall introduce corrected meson fields which minimize the corrected 
energy (||) instead of the mean- field energy (JTJ) . 



3 Corrected meson profiles 

In sect. 2 we considered the static equation of motion (|) and its solution 6(r), in 
the following called "uncorrected" or "static" profile. The energy corrections ([10| - 
[T2| ) were calculated afterwards using the uncorrected profiles. This approach cor- 
responds to a projection after variation. In the present section we want to perform 
a variation after projection by minimizing the corrected energy functional (01) 



5S(r 



-^corr [ 



= (14) 



from which we obtain "corrected" profiles ®J orr {r). The corrected profiles deviate 
from the static ones because the inertial parameters M mert , / and the expectation 
value (-P 2 ) hh depend on the meson profile. Minimizing the corrected energy func- 
tional one allows them to acquire a value which corresponds to a lower corrected 
energy (self-consistent pushing and cranking). Since E^ orr depends explicitly on 
the isospin quantum number T the corrected profiles are different for nucleons 
and A isobars. 

The lower part of fig. 2 illustrates the difference between the self-consistently 
determined moments of inertia for nucleons and A isobars and compares them 
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with the moment of the static hedgehog soliton. Above M ~ 380 MeV, the mo- 
ment for the A isobar is roughly 15% larger than for the nucleon. The harmonic 
average between both momenta agrees quite well the static moment (cf. sect. 4). 
At smaller constituent quark masses, the mean field has a shallow form and the 
valence quarks rather weakly bound. In this case the valence quarks are allowed 
to travel to large distances from the center of the soliton where they are strongly 
affected by the isospin dependent rotational terms in the corrected effective en- 
ergy ([|). The resulting moments of inertia are large and remarkably different 
for nucleon and A isobar. The root-mean square (r. m. s.) radius shows a similar 
behavior (fig. 7). At M < 350 MeV, the uncorrected mean field is too shallow to 
produce a bound quark state and a self-consistent solitonic solution of the static 
equation of motion (|6]) does not exist. After subtracting the translational zero- 
point energy (10) the soliton stabilizes and a self-consistent field configuration 
exists up to M^300 MeV. 
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Fig. 2: Lower part: 
Self-consistently determined iso- 
rotational moments of inertia / 
(full lines) of nucleon and A iso- 
bar in comparison to the static 
moment of the hedgehog soli- 
ton (broken line) as a function 
of the constituent quark mass 
M. The triangles indicate the 
harmonic average of the self- 
consistent moments of both par- 
ticles. 

Upper part: Lagrange parameter 
according to eq. (|T^ ) for nu- 
cleon and A isobar as a function 
of the constituent quark mass. 
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The upper part of fig. 2 shows the Lagrange parameter Q 2 , which is, according 
to eq. (|T3|), necessary to give the soliton the correct isospin. Since the moment of 
inertia decreases with increasing constituent quark mass an increasing rotational 
frequency is needed to reproduce the same isospin. Apart from some minor 
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deviations the Lagrange parameters for nucleons and A isobars differ only in the 
sign. That is why the expectation value (^ 2 ) hh li es between the corresponding 
values for nucleon and isobar. The slight deviations result from the different 
moments of inertia for either particle. 



The corrected equation of motion which follows from condition (II) is much 



more involved than eq. (§) since the correction terms ([ToHT2|) depend on the 



meson profile. Therefore, we apply a variation method in order to find a profile 
function Qj orr (r) which minimizes the corrected soliton energy @. For that 
aim we represent the profile function by a set of numbers 0j, (i = 1, . . . , N) and 
rh n . The 0« = 0(r*) are the values of O(r) at N appropriately selected points 
ri, . . . ,rjy (knots) and m n characterizes the behavior at r^oo. Between the ri 
and r/v the function 0(r) will be reproduced by means of a spline interpolation 
using the 0j's. The asymptotic behavior at small (r<r x ) and large radii (r>rjy) 
requires particular care. Expanding the profile function of the static soliton in 
a power series one can show that the quadratic term vanishes at small r. The 
same is true for the corrected profile function where only the slope is modified by 
the translational zero-point correction. Rotational corrections do not affect the 
behavior at small radii. Therefore we use a linear ansatz at r < r\ in agreement 
with the boundary conditions (|7|) 

( r ) ^ _ n + (Q 1 + n )L. (15) 

n 

The slope is determined by the variation parameter X . 

The behavior at large separations can be determined analytically |17, 



Applying the gradient or heat-kernel expansion on the effective action of the chiral 
soliton in mean-field approximation one can show that the dominating terms agree 
with the non-linear chiral a model with an additional centrifugal force. The 
centrifugal term modifies the asymptotic behavior of the meson profile, which 
is proportional to e~ m,rr at large r, in dependence on the rotational frequency 
Q. In the plane perpendicular to the rotational axis, the pion mass m n , which 
characterizes the exponential descent of the meson field, has to be replaced by a 
modified pion mass pO, [21], p|J 



m w = \jm\ — fi 2 . (16) 

Since Q 2 is negative for the nucleon and positive for the A isobar (cf. fig. 2) 
the nucleon is somewhat slimmer (prolate) than the static hedgehog while the 
isobar is a bit fatter (oblate). Along the cranking axis the soliton is not affected 
by the rotation. Considering this effect in first order we neglect the deviation 
from spherical symmetry and introduce a common effective descent parameter 
fh v (effective pion mass) for all directions which is considered as an additional 
variation parameter. The behavior at r > is then determined by 

r N \ 2 



s 



Restricting ourselves to spherically symmetric solitons we average the effect of 
the rotation over all directions of the rotational axis. Since only two of the three 
spatial directions are affected by the rotation we expect roughly m 2 « m 2 — | f2 2 
instead of eq. (|16|) . Rotational frequencies Q 2 > m 2 lead to negative values of m 2 
in eq. (|T6|). The corresponding 7r field oscillates perpendicular to the cranking 
axis and describes the emission of pions by the soliton. The basis we use for the 
description of wave functions is defined within a finite box [2v| and hence not 
appropriate for the description of oscillating tt fields. The effective pion mass 
resulting from the numerical procedure is close to zero, but real in all cases. In 
this way the basis prevents the soliton from emitting pions. 
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Fig. 3: Square of the effec- 
tive pion mass m n , which 
characterizes the asymptotics 
of the meson field accord- 
ing to eq. (|i7D, as obtained 
from minimizing (||), 
as a function of O 2 . Dif- 
ferent values f2 2 correspond 
to different constituent quark 
masses according to fig. 2. 
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fit performed in the region 
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Fig. 3 displays the numerically determined effective pion masses as a function 
of Q 2 . A fit in the region —6m 2 < fi 2 < m 2 yields m 2 « 1.14m 2 — 0.56 tt 2 in 
fair agreement with our suggestion above. Interpolating to Q 2 = we get an 
effective pion mass which deviates from the pion rest mass by 7 percent. A similar 
deviation (5-10 percent) from the pion rest mass was obtained in a numerical 
variation without any rotational corrections (only translational corrections or no 
corrections at all). This demonstrates the accuracy of our procedure. The obvious 
deviation from the straight line at larger positive values of Q 2 in fig. 3 reveals the 
limitation to descending pion field in our numerical calculation. 

Representing the profile function 0(r) by the set of parameters Oj, (i = 
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1, . . . , N) and rh n the effective energies ([I]) and (|9|) are ordinary functions of 
N +1 variables which can be minimized by standard methods. We tested the 
method by minimizing the static mean-field energy where the self-consistent pro- 
file 0(r) can alternatively be determined in the traditional manner by solving the 
equation of motion (^) iteratively. Here we found that 7 knots are sufficient for 
an accurate reproduction of the meson profile. Moreover, the optimal position of 



the knots was determined. Details can be found in ref. ||24| . Since the corrected 
profiles do not substantially deviate from the uncorrected ones we used the same 
knots for the representation of 0^ rr (V). 



Fig. 4: Profile functions 
6(r) for the static hedge- 
hog soliton (broken line) in 
comparison with the cor- 
rected profiles for nucleon 
and A isobar (full lines) 
calculated for the consti- 
tuent quark mass M = 
600 MeV. Additionally the 
self-consistent profile for 
a translationally corrected 
soliton without rotational 
correction is shown 
(dashed- dotted line). 
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Fig. 4 illustrates the general features of the modification in the profile func- 
tion caused by the correction terms. We selected the relatively large constituent 
quark mass of 600 MeV in order to get a pronounced effect. Translational cor- 
rections modify the behavior of the meson profile at small radii while rotational 
corrections affect the asymptotics at large radii. The first ones increase the slope 
of the internal linear part and make the soliton smaller. The size of the mean 
field and of the corresponding quark distribution is a result of the balance be- 
tween the diverging Fermi motion of the quarks and the attraction between them. 
Subtracting the translational zero-point energy the balance is disturbed since the 




10 



Fermi motion is reduced and the attraction dominates. The soliton shrinks until 
the increasing Fermi motion balances the attraction. 

Rotational corrections are different in sign for nucleon and A isobar and man- 
ifest themselves at large r. They act like a centrifugal force, which is negative for 
the nucleon and positive for the A isobar. In accordance with eq. ( |T6"D and fig. 2 
the exponential descent is larger for nucleons and smaller for A isobars. Fig. 4 
shows that rotational corrections modifies the meson field already at intermediate 
distances from the center of the soliton (r > 0.6 fm). 
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Fig. 5: Upper left figure: Profile function 0(r) for the static soliton in the mass region 
350 MeV < M < 600 MeV; Upper right figure: Deviation AO(r) of the translation- 
ally corrected profile from the static one; Lower figures: Deviations A6 N,A (r) of the 
translationally and rotationally corrected meson profiles for nucleon and A isobar, re- 
spectively. 

Fig. 5 displays the self-consistent meson profiles and the deviations of the 
corrected profiles from the static ones in the whole mass region. As already 
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shown in ref. || the shape of the static profile is practically independent of the 
constituent quark mass (upper left figure). The distortion of the static meson 
profile by the correction terms turns out to be more significant than its variation 
in dependence on M. The effect of the corrections grows in dependence on the 
constituent quark mass M. Heavier constituent quarks are stronger bound and 
produce smaller solitons. Their center-of-mass energy is larger, while the moment 
of inertia is smaller. Both factors increase the correction term and their influence 
on the meson profile. 

4 Expectation values of the soliton 

calculated with corrected meson profiles 

It is the aim of this paper to evaluate the influence of translational and rotational 
corrections of the self-consistent meson fields and to estimate their effect on ex- 
pectation values of solitonic observables. For that reason we compare expectation 
values of observables which are calculated with corrected and uncorrected me- 
son profiles, respectively. In the present paper we will not consider pushing and 
cranking correction to the observables themselves. A full treatment of symmetry 
restoration would have to include these corrections as well. In several cases they 
have turned out to be rather important [p5| , |26| and therefore must be taken 
care of before the results are compared with experimental data. This will be the 
subject of a forthcoming analysis. 

Generally it can be said that the effect coming from the difference between 
uncorrected and corrected meson profiles grows with increasing constituent quark 
mass M in accordance with the increasing size of the correction terms (cf. fig. 1). 
At small constituent quark masses M the valence quarks give the dominating 
contribution to the expectation value in most cases. Increasing the quark mass 
the Dirac sea gets more and more polarized and its contribution gets larger. 
Covering the mass region 300 MeV < M < 600 MeV we are able to control the 
ratio between the contributions resulting from valence and sea quarks. 

Fig. 6 illustrates the energy of the soliton and its components. The total 
corrected energies for nucleons and A isobars are shown in the upper part. The 
experimental nucleon mass is reached at M« 350 MeV, but the predicted mass 
of the isobar is remarkably smaller than the experimental value in the whole 
mass region. Pushing and cranking corrections in the meson field decrease the 
mass furthermore. The difference between full and broken lines indicates the 
gain of energy which is obtained if one minimizes the corrected soliton energy 
instead of its static value. Here we compare the corrected total soliton energy @ 
calculated for the corrected profiles 0^ rr with the same energy calculated for the 
static profile G. While the total energies do not differ very much from each other 
valence and sea quark energies are rather different in both cases (central and lower 
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part of fig. 6). For the static soliton we have a dominating valence contribution at 
small constituent quark masses (M < 420 MeV). In the energy of the corrected 
soliton, the valence quarks dominate up to larger masses (M < 500 MeV). The 
correction terms even prevent the valence quark energy from getting negative 
what happens at M« 750 MeV in the uncorrected case || . The field energy (fj) 
of the mesons is strongly affected by the correction (up to a factor of 2). However, 
its contribution to the total soliton energy is marginal due to the restriction to 
the chiral circle. 



1 

> 1 

0) 



1 



200 
000 
800 
600 
400 

000 
800 
600 
400 
200 


800 
600 
400 
200 





sea - 



mes 



"val ^^^T^- 


A isobar 1 
^_ . 


sea 


"~ * * m - — . 




300 400 500 

M (MeV) 



600 



Fig. 6: Upper part: To- 
tal corrected energies £^J,rr 
for nucleon and A isobar 
calculated for the static 
meson profile G (broken 
lines) in comparison to the 
same energy calculated for 
the corrected profiles 
(full lines). 

Central and lower part: 
Contributions of the va- 
lence (val) and sea quarks 
(sea), and of the mesonic 
field (mes) to the total 
soliton energy E m { calcu- 
lated for the static meson 
profile (broken lines) and 
for the corrected profiles 
(full lines) as a function of 
the constituent quark mass 
M. 



The effect of the corrections on an expectation value depends crucially on 
the relation between valence and sea-quark contributions. If the valence quarks, 
which are usually well bound inside the soliton, give the dominating contribution, 
the shrinking of the meson profile in the internal region (see fig. 4) is the most 
important effect. It is caused by the center-of-mass correction and is independent 
of the isospin quantum number. Rotational corrections have almost no effect 
on the valence quarks. Sea quarks are either less bound or unbound and hence 
experience the action of rotational corrections which are opposite for nucleons and 
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A isobars. Considering the sea contribution of an expectation value in dependence 
on the radius r one notices that contributions from the internal part of the soliton 
may cancel out contributions from external layers. This is not the case for the 
valence part since the corresponding wave function has no node. Modifications 
in the meson profile influence the complicated interplay between positive and 
negative contributions and may affect the total sea contribution considerably. As 
a result sign and size of the effect depend strongly on the observable. In the 
following we give some examples in order to illustrate what might happen. 
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Fig. 7: A-nucleon mass 
splitting A£an , isovec- 
tor magnetic moment 
Li T=1 = fi p — axial- 
vector coupling strength 
<7A and isoscalar root 
mean-squared (r. m. s.) 
baryon radius R calcu- 
lated for the translation- 
ally and rotationally cor- 
rected nucleon (N) and 
A profiles (full lines) in 
comparison to the same 
quantities for the static 
soliton (broken lines), as 
a function of the con- 
stituent quark mass M. 



The (isoscalar) r. m. s. radius R of the quark distribution (lower part of fig. 7) 
is dominated by the valence quarks by nearly 95 percent. Hence nucleon and A 
isobar have almost the same radius. The reduction of the radius with respect 
to the static hedgehog by a few percent is in quantitative agreement with the 
steeper behavior of the corrected meson profile in the central region of the soliton 
(cf. figs. 4,5). For very small constituent quark masses the valence quarks are so 
weakly bound that they may reach larger separations from the center. Here they 
are exposed to the rotational corrections which lead to a considerable shrinking for 
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the nucleon and a remarkable swelling for the A isobar with respect to the static 
hedgehog. The r. m. s. radius is particularly sensitive to the quark distribution at 
large radii since the quark density is weighted by a factor r 4 . 

Now let us investigate the axial-vector coupling constant gx of the nucleon 
(fig. 7). In order to study the pure influence of the corrected meson field we 
consider the expectation value of the same operator || without any rotational 
corrections. This uncorrected observable accounts for roughly 65 percent of the 
experimental value of g& [25, ^6|. Apart from rotational corrections the coupling 
constant for the A isobar differs from the nucleon by an additional factor which 
results from the different spin-isospin structure of both particles. The figure shows 
the #a of the proton calculated for meson profiles which have been corrected for 
both nucleons and A isobars. Here the dominance of the valence quarks is less 
pronounced as in the case of R and the different asymptotic behavior of nucleon 
and A profiles results in slightly different expectation values. Roughly spoken 
translational and rotational corrections cancel each other for the nucleon while 
they add up (negatively) for the A isobar. 

In the isovector magnetic moment of the nucleon (fi^ =1 ={i p — /i n ), contribu- 
tions from larger radii are favored by a factor r 3 |l], |26|. Sea contributions amount 
to 20. . . 30 percent of the total value. This leads to a noticeable difference be- 
tween the expectation values for nucleon and A profiles. Here, the correction is 
negative for the nucleon and almost cancels out for the A isobar. 

In the last example we consider the delta- nucleon mass splitting AE^n- In 
the case of the static soliton the splitting results solely from the different cranking 
energies of both particles and is uniquely determined by the moment of inertia / 
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Using the corrected meson profiles, which are slightly different for nucleons and 
A isobars, the self-consistent values of static energy (H), translational zero-point 
energy ( |T0| ) and of the moment of inertia are different as well. Now the mass 
difference is given by 
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where the additional upper index N, A indicates the corresponding particle. In 
the upper part of fig. 7 we see that the simple expression (|18|) with the static 
moment of inertia agrees with the more complicated expression (fl9|) for corrected 
profiles very well. The analysis of the various terms in eq. (|19"D shows that the 
first two differences are rather small (~ lOMeV). In spite of slightly different 
moments (cf. fig. 2) of nucleon and A isobar the last term is nearly identical with 
the result of eq. (113) . That means that the harmonic average defined by 
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agrees with the moment of the static hedgehog soliton (cf. fig. 2). The mass 
splitting provides an example of an almost complete canceling of meson-field 
corrections. 

The magnitude of the evaluated effects does not exceed a level which de- 
scribes the reliability of the model in particular in the physically relevant region 
of small constituent quark masses. The situation may be different if one consid- 
ers observables with oscillating radial dependence, e.g. formfactors. Here a minor 
modification in the meson profile may have a larger effect. 

5 Conclusions 

We have evaluated pushing and cranking corrections on solitonic meson fields 
of the bosonized NJL model in one-loop approximation and studied their effect 
on the quark distribution and other observables of the soliton. The corrected 
fields have been obtained by minimizing an energy functional which differs from 
the mean-field energy by several correction terms. These terms have been in- 
troduced to remove the energy of spurious translational and rotational modes 
and to equip the soliton with the quantum numbers of nucleon or A isobar on 
a semi-classical level. The meson fields are restricted to the chiral circle and to 
hedgehog configurations with a modified asymptotic behavior. 

We have studied solitons and their expectation values in the region 300 MeV < 
M < 600 MeV of the constituent quark masses M. The results illustrate the re- 
sponse of the meson field to the corrections and quantify their effect on expecta- 
tion values. Despite the big correction terms (up to 50 percent of the total soliton 
energy) meson and quark fields are only moderately affected. The (approximate) 
restoration of translational symmetry increases the slope of the meson profile in 
the central region and makes the soliton smaller. This effect prevents the valence 
quarks from diving into the negative-energy region up to very large constituent 
quark masses. It stabilizes the soliton by reducing the kinetic energy of the 
quarks. As a result, a stable soliton exists at quark masses below 350MeV(up to 
300 MeV) while the uncorrected soliton is heavier than 3 free quarks and de- 
cays, in this mass region. Rotational corrections affect the asymptotic behavior of 
the fields at large radii. They depend on the spin and isospin quantum numbers 
and therefore have different effects on nucleon and A isobar. Both translational 
and rotational corrections grow with increasing constituent quark mass. 

Expectation values of the soliton are affected by the corrections to various ex- 
tents. The isoscalar r. m. s. radius of the nucleon is reduced by only a few percent 
(~ 3-5 percent). Larger changes can be observed if one considers valence and 
sea quark contributions separately. Axial-vector coupling strength and isovector 
magnetic moment are modified up to 40 percent. In the physically relevant region 
of small constituent quark masses (350 MeV < M < 450 MeV), the corrections do 
not exceed 20 percent. The A- nucleon mass splitting is practically not affected 
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by the corrections despite the slightly different moments of inertia calculated for 
nucleon and A isobar. 

The experimental values of r. m. s. radius and AN mass splitting are well re- 
produced for light constituent quark masses around 400 MeV. Here, the effect 
of translational and rotational corrections to the mean field is small and does 
not exceed the level of accuracy of the model and its numerical treatment. The 
quantities g^ and yU T=1 are stronger affected by the corrections. Up to now the 
theoretical values underestimate the experimental ones significantly. The calcu- 
lated values are even smaller with the corrections than without. However it must 
be noted that our treatment is incomplete, because in order to be able to compare 
with the experiment one has to take into account that there are also rotational 
corrections to the corresponding operators, which have turned out be quite large 
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